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Cosmic balloons are spherical domain walls with relativistic particles trapped 


inside. We derive the exact mass and radius relations for a static cosmic bal- 
loon using Gauss- Codazzi equations. The cosmic balloon mass as a function 
of its radius, M(R), is found to have a functional form similar to that of 
fermion soliton stars, with a fixed point at 2 GM(R)/R ~ 0.486, which corre- 
sponds to the limit of infinite central density. We derive a simple analytical 
approximation for the mass density of a spherically symmetric relativistic gas 
star. When applied to the computation of the mass and radius of a cosmic 
balloon, the analytical approximation yields fairly good agreement with the 
exact numerical solutions. 
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1. Introduction 


A particle with a discrete symmetry may have a different mass on either side of the 
domain wall, which may form in an early Universe phase transition with the spontaneous 
breaking of the discrete symmetry. When a particle has less energy than the rest mass 
it would have on the other side of the wall, it can become trapped inside a region with a 
closed domain wall as the boundary. Closed domain walls filled with such particles have 
been named cosmic balloons” [1]. Due to effects such as the emission of gravitational 
radiation, irregular shapes of balloons oscillate and eventually settle down to spheres. 
Here we are only concerned with a spherical cosmic balloon. 

Cosmic balloons are stable objects. The pressure of the relativistic gas inside the 
balloon balances the surface tension of the wall. An incarnation of a cosmic balloon is 
a neutrino ball [2], which is a spherical domain wall with light right-handed neutrinos 
trapped inside. Details of neutrino balls have been worked out which seem to indicate that 
they are plausible astrophysical objects [2], possibly providing an alternative explaination 
for the mass of quasars and other astrophysical phenomena. 

In this work, we follow Ref[l] in the study of general properties of cosmic balloons, 
in particular, the mass of the balloon as a function of its radius. The analytical ap- 
proximation in Ref[l] has an extremely involved numerical form, and the mass versus 
radius curve found there is numerically incomplete. In another previous work with sim- 
ilar objectives, Ref[3], an analytical approximation was applied beyond its valid range 
and led to incorrect numerical solutions of the mass and radius. In this work, we find 
simple and transparent analytical expressions which are useful in helping us understand 
the cosmic balloon solutions qualitatively, and we compute the complete set of mass and 
radius numerically with and without the application of the analytical approximation. 
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2. Mass and radius relations 


To find the radius R and mass M of a cosmic balloon, we study the static configuration 

of a spherical domain wall containing gas. Let the metric inside the domain wall (r < R) 
be 

dr 2 = B(r)dt 2 - A(r) dr 2 - r 2 {d9 2 + sin 2 0# 2 ). (2 


The motion of spherical domain walls containing vacuum has been solved in Ref.(4), 
using the Gauss-Codazzi formalism. Following Ref.(4), we find the equations of motion 
for a spherical domain wall containing a perfect fluid with pressure p and density p: 


(a + /3)R + 


aGM 

R 2 



(a — P)R + 


aGM 

R 2 


£ 

2 


A ' R 2 4 - B ' 

-R +-a 

T f>2 . B' 2 

A R + ~B a 


—a/3 


2(a + /?) 


R 


2p 2 R 2 

A(p + p ) 


—4a/37rGcr, 


( 2 . 2 ) 


Here a dot denotes a derivative with respect to proper time, and a prime denotes a 
derivative with respect to r — R. We have defined 


as [/T'C^ + R 2 ] 172 , 
a r. 2 GM . 2 i‘/ 2 

^l 1 - — + iJ ] 


(2.3) 


Af is the total mass of a cosmic balloon of radius R. 

The equation for hydrostatic equilibrium is [5] 

B[ _ 2 p f 

B p + p' (2.4) 

This and the Einstein equations can be combined to give a single differentia] equation: 
- r 2 p'(r) = GM{r)p{r) 


fi + PWl 

1 t 47r r 3 p(r) 

' 2 GM(r)' 

P{r)\ 

M(r) 

1 — — 

r 


-1 


(2.5) 


where we have defined M(r) such that 


M'(r) = 4ir r 2 p(r), M(0) = 0. 


( 2 . 6 ) 
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For 4(0) finite, 


A(r) = 


1 - 


2 GM(r) 


i -l 


(2.7) 


Using the above well known results in Eqs.(2.2) and setting R = 0 = R, we find 

M{R) = M(R ) + 2irR 2 a(a + /?), ^ 

(2.9) 




with 


a = 


2 2 GM(R) 

R 


1 1/2 


/? = 


1 - 


2GM1 1/2 




( 2 . 10 ) 


Eq.(2.8) can be rewritten to give 


0 = ot — iirRGa, for R < — / oin 

AnGa 

A cosmic balloon with R > a/(4? xGa) is contained within its Schwarzschild radius, a 
case in which we are not interested here. Hence 


M{R) — M(R) + 4t rR 2 a[a - 2tt RGa]. (2 12) 

Note that M(R) is the total energy of the matter and the gravitational field inside the 
cosmic balloon [5]. Eq.(2.9) gives the radius R of a cosmic balloon as a function of its 
surface tension cr, it can be rewritten as 

/ m 2a<r GaM(R) a 0 

p{R) - — + ' ^ + j 3 2 ttGo 2 . (2.13) 

We see that the gas pressure at R is balanced by the surface tension of the domain wall, 

the gravitational force on the wall due to the total mass inside the wall and the wall 
itself. 
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3. An exact solution 


Let us define dimensionless variables as follows: 


x(r) 

f 

p*(r) 

M(R) 


2 GM(r) 


V 

p(r) 

Po 

M(R) 

M 0 


( r o S [ 487 rG<r] -1 ) , 


(/>o = [8 tt G rj*]- 1 ) , 

(*■ 3 )- 


Note that the length scale r 0 
ro ~ 1.95 x 10 I3 cm 

M 0 6.6 x 1O 7 M 0 


and the mass scale Mo can be written as 

a I ' 1 
(TeV) 3 J ’ 

<7 I" 1 

(TeV) 3 _ ' 


(3.1) 


(3.2) 


This is interesting since both R and M(R) have maxima on the order of 1, as we shall 
see. a 1 / 3 is usually with the sca , e of the discrete symmetry breaking) but jt 

could be made anything depending on the particular model of cosmic balloon. 

Eq.(2.5) can be rewritten as 


= -f rh t *' 2 + 5x ' fi ^ + - vm?)} - 


where the primes denote differentiation with respect to r, and 

H 1 * 1 ^ 

III 

< 

/i(0) = 0 , 

Mr) = £ 
r z 

MO) = , - 

3 

3 

II 
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II 

II 

O 

if p( 0 ) is finite. 


(3.3) 


(3.4) 


Note that p( 0) dependence only comes in through the boundary condition of the function 
/ 2 (r). In terms of x and r, 


P*( r ) = Z 2 [% + r x '] . 


(3.5) 
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The scaled radius of the cosmic balloon R is given by the intersection of p*(r) and 

P ' (R) = 35 ( 3a + ■ (3.6) 

The scaled mass of the cosmic balloon is 


M(R) = R 


,-kv , R ( R 

x{R) + — a 

6 V 24 


(3.7) 


Here 


a = [l -x{R)] l/2 , 


P = a- R/ 12, 


for R < 12a. 


(3.8) 


Eq.(3.3) has an obvious solution [5] 


3 

x = - 

V 


P ^ 56ttGV 2 ° r p ’^ ~ 77=2’ 


which has infinite central mass density. The corresponding cosmic balloon radius and 


mass are: 


^ 4 


R = (5 - V22) ~ 0.468, 

M(R) = ^(^2 ~ 25-y /22) 


63 V7 


0.2275, 


2 GM(R) _ M(R) _ 2(20 - y/22) 
R R ~ 63 

4. Numerical solutions 


~ 0.486. 


(3.9) 


The numerical solutions of the cosmic balloon mass function can be obtained in a straight- 
forward way. It is most convenient to use the parameter q to track the solutions of 
Eq.(3.3). Note that 

_ SnGp(0)r^ 


1 = ^^ocp(0), 


i>*(0) = 3 q. 


(4.1) 


For a given value of q, Eq.(3.3) can be integrated to find x(r) and hence p*(r). p*(r) and 
p*(R) = (3a + p~ l )/4 R have two intersections for q > q c , which give two values of the 
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cosmic balloon radius R , corresponding to two different branches of the cosmic balloon 
mass/radius solutions (see Fig. 1). The two branches of solutions meet at q = q c ~ 1, 
for q < q c (or p( 0) < p c (0)) there are no solutions and cosmic balloons do not exist. q c 


corresponds to the minimum of the central mass density for a cosmic balloon. 

Fig.2 and Fig.3 are the radius R and mass M(R) versus q respectively. For increasing 
q (i.e., p(0)), the first branch of solutions gives decreasing R and M(R); the second branch 
of solutions gives R and M{R) with complex and similar behavior. In the second branch 
of solutions, R increases until it reaches the maximum ( R m ~ 0.732), then it decreases 
until it reaches the minimum, it then increases again until it reaches a second (smaller) 
maximum, then it decreases again to a second (larger) minimum, and continues in this 
manner. From Fig. 2, it is clear that R oscillates with damping amplitude around 
R ~ 0.468 for large q . M(R) has the same behavior but with a phase lag in ?[p(0)] 
compared to R ( M m cx 0.374, R^ m ~ 0.71), and it oscillates around ~M(R) ~ 0.2275 
for large q (see Fig. 3). This is not surprising, since R ~ 0.468 and M(R) ~ 0.2275 
correspond to the exact solution of p(r) in the limit of p(0) -► oo (i.e., q oo). In the 
plot of M(R) versus R (Fig. 4), we see that the mass function spirals in toward the fixed 
point (R ~ 0.468, M(R) ~ 0.2275), which satisfies 2GM(R)/R ~ 0.486. 

The maximum mass and the corresponding radius for a cosmic balloon are 


M max = M m M 0 ~ 2.47 x 10 7 M. 


<7 


0 


(TeV) 3 J 




R-Mmax — Rmtu ^0 — 1.38 X 1 0 ^ 'Vlll 

(TeV) 3 J ’ 

Since 2 GM(R)/R = M(R)/R , the fractional red shift ^ of a spectral line emitted 
the surface of a cosmic balloon [5] 


(4.2) 

from 



/ 2 GM\-V 2 

l 1 -- T ) ~ x 


(4.3) 


is independent of the actual physical scales r 0 and M 0 , i.e., independent of a (the sur- 
face tension of the domain wall). The exact numerical solution shows that 2GM/R 
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reaches maximum ((2 GM/R) max ~ 0.56) after the mass has already reached maxi- 
mum. Before M = M max , 2GM/R increases with increasing M. The maximum of 
mass then corresponds to the maximum of the fractional red shift 2 for a stable cosmic 

balloon, z max ~ 0.46 (2 GM m /R Mm ~ 0.53). At the minimum of the central mass density 
(R c ~ 0.56, M c ~ 0.17, 2 GM C /R C ~ 0.30), 2 C ~ 0.20. 


5. Analytical approximations 

To obtain useful analytical approximations, we scale Eq.(3.3) such that its solutions are 
independent of the central mass density p( 0) (i.e., independent of q). Define 

t = y^Sgr. 


Eq.(3.3) becomes 

f + 5if w + ( 7x - 3 )/ 2 *m] , 

where the dots denote differentiation with respect to t, and 


m = f . m = 0, 

m = fm = i 

x(t ~ 0) = 0 = x(t = 0). 


Note that Eq.(5.1) is exactly the same in form as Eq.(3.3), but without p( 0) dependence 
in the boundary conditions. Correspondingly, 




m = h x+ti ]• 


(5.3) 


Combining Eqs.(5.3) and (5.1), we can express x(t) in terms of p(t): 


[(p/p) +2pt/3]t 


CM 

1 

t-i 

(5.4) 
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The radius of the cosmic balloon is now given by 

7{tR) = vfe ( 3o(<r) + ]9(b) ’ (5.5) 

where t R = y/ 3 qR, and 


a(t R ) = [1 -x(t R )] 1/2 , 


P(tn) = <*(*/?) - 


tR 


I2y/3q' (5.6) 

Given the function p(t), we can find R(q). The mass of a cosmic balloon is (see Eq.(3.7)) 

r 15 / "n \ 1 

M(R) = R 


+ 1 U(t R ) - ~ 


(5.7) 


Recall that M(R) and R are defined in Eq.(3.1). 
The analytical solution from Ref[3] 


PiW = 


cosh 2 (A;£) ’ 


2 

k=di 


(5.8) 


is an extremely good approximation for t < 1. The corresponding x(t ) is 

/ _ (kt)tanh(kt) - t 2 /[ 3 cosh 2 (kt )] 

(kt) tanh(fct) + 1 ’ (5-9) 

Eq.(5.1) still has the exact solution x = 3/7, which has the wrong boundary conditions 
at ^ = 0. However, the solution with the correct boundary conditions at t = 0 should 
approach x = 3/7 in the limit of large t, as indicated by our exact numerical solutions. 
To the second order of a small parameter A, the perturbation around x = 3/7 gives 

*2(«) = f + Jjl COSp(t)] + 

7 A 2 f 139 19 _ x 

32 \ -7 + 54- cos[20(f)J - — V47 sin[2#(i)] j , (5.10) 

where 


e(t) = 




+ tan 1 



7T 

B' 


(5.11) 
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A and B are constants which are independent of p( 0). A good fit of Eq.(5.10) with the 
exact numerical solution gives 


A ~ —0.18, 


B ~ 3.54. 


(5.12) 


Eq.(5.3) gives the corresponding p(t) for x 2 (t): 


p 2 (t) = 


t 2 


7 + 
7 A 2 


~ (cos[0(f)]-^sin[0(t)]) 


— cos[2i9(f )] — — \/47 sin 


(5.13) 


32 ( 3/2 \ 2 

Since Eq.(5.8) is valid for small /, while Eq.(5.13) is valid for large t , we can construct an 
approximate solution by combining the two solutions with appropriate weight functions, 


i.e. 


p{t) — Wi(t) pi(t) + w 2 (t) p 2 (t), 
where Wi(£) and w 2 (t) are weight functions. For example 


Wi(t) = 


1 + e f / ni t n 2 ’ 
A good choice is 


W 2 (t) = 


t n 2 




wi = 20, 


n 2 = 8. 


(5.14) 


(5.15) 


(5.16) 


The cosmic balloon radius as a function of the central mass density p(0) can be 
computed by using Eqs.(5.14) and (5.5). We can write 

x(t) = w,(t) Xl (t) + w 2 (t)x 2 (t). 

The ma^s of the cosmic balloon is given by Eq.(5.7). 

In Fig.l, we see that Eq.(5.14) is an extremely good approximation of the exact mass 
density p(t) (scaled to be independent of p(0)). The radius R and mass M(R) obtained 

by using Eqs.(5.14) and (5.17) are plotted in dashed lines in Figs. 2-4, they are quite 
close to the exact solutions. 
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6. Remarks 


The spiral behaviour in M(R) versus R of cosmic balloons (see Fig.4) resembles that of 
fermion stars [6], This is not surprising since cosmic balloons and soliton stars are similar 
objects. dM/dp( 0) = 0 generally signifies the change from stability to instability (or vice 
versa), and dM/dp( 0) > 0 indicates stability. However, due to dynamical complications, 
the range of M(R) and R beyond the first mass maximum (dM/dp( 0) = 0) on the curve 
in Fig. 4 does not correspond to stable configurations of cosmic balloons, although the 
conditions for stability are satisfied each time the curve bends upward (for increasing 
p(0)). It is possible that the existence of the fixed point (R ~ 0.468, M(R) ~ 0.2275 
2 GM{R)/R ~ 0.486) indicates that the cosmic balloons beyond the stability range os- 
cillate around a finite configuration with 2 GM(R)/R ~ 0.486, instead of collapsing into 
black holes. 
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Figure Captions 


Fig. 1. p{t) and p(t R ) = [3a(f K ) + p-\t R )]/ (4^ qt R ), with t = ^ and q = 
87rGp(0)^/3 = 0.9,2,10 (dotted lines from right to left). The dashed line is the analytical 

approximation of p(t) from Eq.(5.14), it almost overlaps with the exact solution (the solid 
line). 

Fig. 2. R versus q = 87rGp(0)r^/3. The solid line is the exact solution, while the 
dashed line is the result of analytical approximation. The dotted line is R = 0.468. 

Fig. 3. M(R) versus q = 8xGp(0)rg/3. The solid line is the exact solution, while the 
dashed line is the result of analytical approximation. The dotted line is M(R) = 0.2275. 

Fig. 4. M{R) versus R. The solid line is the exact solution, while the dashed line is 
the result of analytical approximation. 
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